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ABSTRACT 

Vertical  wave  barriers,  also  known  as  vertical  wall  breakwaters,  are  designed  to 
reduce  the  amount  of  wave  energy  behind  the  barrier  to  allow  the  safe  mooring  of  vessels. 
These  structures  are  constructed  of  wooden  or  concrete  panels  mounted  to  piles.  The  panels 
extend  from  above  the  water  surface  down  to  approximately  mid-depth,  leaving  a  gap  near 
the  sea  floor.  The  existing  procedures  used  to  design  vertical  wave  barriers  remain 
essentially  unchanged  since  1960.  These  procedures  produce  uncertainties  when  used  to 
predict  transmitted  wave  height  and  wave  loading  (force)  on  the  barriers;  it  is  hoped  that 
new  approaches  developed  in  this  project  may  be  adopted  into  modem  design  guidelines. 

The  theoretical  approach  of  this  investigation  involved  modification  of  existing  linear 
wave  theory  to  account  for  the  effects  of  reflected  waves.  This  resulted  in  a  new  theory  for 
both  wave  transmission  and  for  wave  forces  on  the  barrier.  Wave  transmission  and  reflection 
and  resultant  forces  were  measured  when  a  model  wave  barrier  was  subjected  to  regular  and 
irregular  (random)  waves.  Additionally,  flow  characteristics  underneath  the  breakwater  were 
measured.  With  the  recorded  data  the  new  theories  can  be  verified  and  substantial  forward 
progress  can  be  made  in  the  design  procedures  for  these  barriers. 
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INTRODUCTION 

Recently,  there  has  emerged  a  renewed  interest  among  ocean  engineers  in  vertical, 
thin  “wave  barriers”  as  a  lower-cost  breakwater  option  to  protect  small  recreation  and 
fishing  harbors.  These  breakwaters  are  generally  constructed  of  sheets  or  panels  of  rigid 
material  attached  to  vertical  piles  which  are  embedded  in  the  ocean  floor,  as  shown  in  Figure 
1.  The  sheets  of  material  extend  only  partially  to  the  sea  floor,  leaving  a  gap  between  the 
bottom  of  the  breakwater  and  the  sea  floor.  Unlike  conventional  rubble  mound  breakwaters, 
which  are  built  with  a  trapezoidal  cross-section,  wave  barriers  or  wave  “screens”  are 
economically  feasible  in  intermediate  depths  of  water  ranging  from  30  to  60  feet.  This  type 
of  breakwater  is  well-suited  for  wave  protection  in  areas  subject  to  waves  of  about  three  feet 
(one  meter)  in  height  with  periods  of  less  than  six  seconds.  A  vertical  wall  breakwater  can 
typically  be  built  for  about  1/3  of  the  $10,000  per  foot  that  rubble  mound  breakwaters  can 
cost.  (Gilman  and  Nottingham,  1992).  A  primary  advantage  of  vertical  wave  barriers  is  that 
water  can  circulate  underneath  the  breakwater,  helping  to  maintain  harbor  cleanliness  and 
water  quality.  This  design  does  not  completely  prevent  wave  transmission,  however,  so  care 
must  be  taken  to  extend  the  wall  to  sufficient  depth  to  limit  wave  heights  within  the  harbor 
to  an  acceptable  level. 

The  specific  design  and  material  composition  of  wave  barriers  depends  on  the 
decisions  of  the  designing  engineer  when  taking  the  local  environment,  application,  and  cost 
into  consideration.  In  the  Seattle  Waterfront  Project  (Gilman,  1995)  for  instance,  the  wave 
barrier  is  constructed  of  concrete  panels  that  are  attached  to  vertical  steel  piles.  These 
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Figure  1.  Diagram  of  vertical  wave  barrier  from  Peratrovich,  Nottingham,  and  Drage,  Inc., 
(Issue  #39,  undated). 
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vertical  piles  are,  in  turn,  connected  to  a  pier  composed  of  vertical  and  batter  piles  and  a 
concrete  cap  as  shown  in  Figure  2.  This  design  occupies  minimal  area  within  the  harbor, 
supports  a  pier  deck,  and  allows  ships  to  be  tied-up  alongside  the  breakwater  [Peratrovich, 
Nottingham,  and  Drage,  Inc.  (PN&D),  Issue  #39,  undated).  The  Coast  Guard’s  Tillamook 
Station  in  Garibaldi,  Oregon,  is  protected  by  a  combination  steel  and  wood  breakwater.  This 
breakwater  is  designed  to  withstand  waves  with  periods  up  to  five  seconds  and  is  actually 
slanted  at  an  angle  to  the  sea  floor.  (Gilman  and  Nottingham,  1 992).  A  wave  barrier  similar 
to  the  Tillamook  barrier  is  shown  in  Figure  3. 

Vertical  wave  barriers  were  conceived  based  on  the  theory  that  most  of  the  energy 
of  a  wave  is  concentrated  near  the  water  surface.  The  kinetic  energy  within  a  wave  rapidly 
decays  as  depth  increases.  The  fluid  particle  velocities,  and  hence  wave  energy,  decrease  to 
nearly  zero  at  a  point  one-half  of  the  wave  length  below  the  surface  of  the  water,  a  depth 
sometimes  referred  to  as  the  “wave  base”  (Figure  4).  If  a  wave  has  a  wave  base  deeper  than 
the  depth  that  the  breakwater  extends  beneath  the  water  surface  (penetration),  then  some 
amount  of  wave  energy  will  pass  under  the  breakwater.  The  amount  of  transmission  depends 
on  the  disparity  between  these  two  values  of  breakwater  penetration  and  one-half  of  the 
wave  length;  the  greater  the  difference,  the  greater  the  transmission. 

There  was  considerable  research  done  on  wave  barriers  in  the  1950s  and  ‘60s. 
However,  they  have  not  been  widely  used  or  studied  in  the  past  30  years.  Nearly  all  existing 
theory  concerning  the  interaction  of  waves  with  partial  vertical  wall  breakwaters  is  derived 
from  a  paper  by  Wiegel  (1960)  which  outlines  a  simplified  theory  for  the  transmission  of 
waves  by  these  breakwaters.  The  theory  to  predict  wave  forces  on  wave  barriers  was 


Figure  3.  Side  view  of  wave  barrier  similar  to  that  at  Tillamook  Coast  Guard  Station, 
PN&D  (Issue  #30,  undated). 


Deep  Water 


Figure  4.  Diagram  of  deep  water  wave  base.  Water  particle  orbits  and  wave  energy 
decrease  with  depth.  Sorenson  (1993),  pp.  16 
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developed  by  Sainflou  (1928)  and  Miche  (1944).  Rundgren  (1958)  made  experimental 
observations  and  later  modified  Miche’s  higher-order  wave  theory.  The  Sainflou  (1928)  and 
Miche-Rundgren  (1958)  theories  are  currently  used  by  the  Army  Corps  of  Engineers  (1984) 
in  their  Shore  Protection  Manual.  This  manual,  and  to  a  lesser  extent  the  Design  Manual 
26.2  of  the  Naval  Facilities  Engineering  Command  (1982),  is  the  standard  for  coastal 
engineering  and  promulgates,  among  other  things,  suggested  design  methods  and  parameters 
for  the  construction  of  vertical  wave  barriers. 

Existing  equations  used  to  determine  the  wave  transmission  past,  and  wave  forces 
on,  wave  barriers  are  believed  to  be  overly  conservative.  Specifically,  it  is  believed  that  the 
method  of  predicting  wave  transmission  as  developed  by  Wiegel  (1960)  over-estimates  wave 
transmission  for  most  conditions  of  interest.  The  Miche-Rundgren  method  for  computing 
force,  as  outlined  in  NAVFAC  DM  26.2,  assumes  the  worst-case  condition  of  a  wave  crest 
on  the  incident  side  of  the  barrier  and  a  wave  trough  on  the  other.  This  maximizes  the 
difference  in  wave  pressure  across  the  barrier  and  is  an  unrealistic  assumption.  The  predicted 
forces  from  this  method  are  greater  than  those  that  would  actually  occur,  which  can  result 
in  the  over  design  of  these  breakwaters.  Larger  cities  such  as  Seattle  can  sometimes  afford 
to  build  these  breakwaters.  However,  the  inflated  costs  of  this  type  of  breakwater  often 
prevent  their  use  in  smaller  communities.  The  Canadian  government,  for  example,  has 
considered  the  use  of  this  type  of  breakwater  to  protect  numerous  small  fishing  harbors  on 
the  British  Columbia  coast,  but  then  usually  rejected  the  design  due  to  excessive  cost  (Wave 
Barrier  Symposium,  1995). 

This  Trident  study  was  undertaken  to  improve  the  theory  used  in  the  design  of  wave 


9 


barriers  and  to  generate  laboratory  data  to  validate  the  new  theory.  This  study  had  four  main 
parts:  (1)  literature  search  of  experiments  and  theories  in  related  areas;  (2)  development  of 
modified  theories  for  wave  forces  on  and  wave  transmission  under  a  thin  wall;  (3)  regular 
and  irregular  wave  tests  on  a  thin,  rigid  vertical  breakwater  to  obtain  experimental  data;  and 
(4)  evaluation  of  theories  based  on  measured  values  of  wave  transmission  and  wave  forces. 

The  literature  search  for  prior  experiments  and  theories  for  transmission  and  forces 
was  conducted  in  order  to  gain  an  idea  of  the  type  of  testing  performed  in  previous 
experiments,  the  theories  that  were  used  and/or  developed  in  these  tests,  real-world 
applications  of  these  breakwaters,  and  the  construction  methods  of  these  breakwaters. 
Modified  theories  concerning  forces  and  wave  transmission  were  developed  based  on 
WiegePs  (1960)  theory  and  differ  mainly  as  a  result  of  the  inclusion  of  wave  reflection  in  the 
development  of  theory.  The  laboratory  data  are  based  on  tests  conducted  in  the  120  foot 
wave/towing  tank  of  the  U.S.  Naval  Academy  Hydromechanics  Laboratory  on  a  model  rigid, 
thin  vertical  wave  barrier.  The  regular  and  irregular  tests  were  conducted  to  measure  wave 
setup  on  the  structure,  the  resultant  forces  on  the  wall  due  to  wave  action,  and  transmitted 
wave  height  and  to  then  compare  the  regular  and  irregular  results  to  the  modified  theories. 
Particle  velocities  under  the  wall  and  pressures  in  front  of  the  wall  were  also  recorded  for 


some  tests. 


The  primary  purpose  of  a  vertical  wave  barrier  is  to  reflect  incident  wave  energy, 
thereby  creating  an  acceptable  wave  climate  in  the  area  behind  the  barrier  acceptable  for  the 
mooring  of  small  boats.  For  most  breakwaters  used  in  harbor-protection,  effectiveness  is 
determined  by  their  reduction  of  transmitted  wave  height.  This  is  quantified  through  the 
coefficient  of  transmission,  K, ,  which  is  equal  to  transmitted  wave  height,  F(  ,  divided  by 
incident  wave  height,  H; . 


The  transmission  allowed  by  wave  barrier  breakwaters  is  variable  and  is  primarily  a  function 
of  wave  length,  water  depth,  and  breakwater  penetration.  Wave  forces,  on  the  other  hand, 
depend  predominantly  on  wave  reflection  and  also  display  a  much  stronger  dependence  on 
wave  height  than  is  shown  by  wave  transmission. 

Wave  period  is  an  important  factor  because  wave  length  is  dependent  solely  on  wave 
period.  In  deep  water,  the  wave  length,  L  (Figure  5),  is  given  by  the  relationship 


(2) 


where  g  is  the  force  of  gravity,  and  T  is  the  wave  period.  It  is  wave  length,  not  wave  height, 
that  has  a  primary  influence  on  transmission.  Seemingly  contrary  to  common  sense,  the 
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amount  of  wave  energy  present  has  no  influence  on  the  coefficient  of  transmission.  The 
vertical  distribution  of  this  energy  is  the  critical  factor,  and  wave  length  is  what  determines 
the  distribution  of  the  energy  in  a  wave.  The  magnitude  of  transmission  is  dependent  on 
how  much  of  this  energy  is  present  at  a  depth  that  will  pass  below  the  vertical  wall.  The 
depth  to  which  water  particles  are  affected  by  a  wave  is  called  the  “wave  base”  and  is  equal 
to  one-half  the  wave  length.  Thus,  the  transmission  coefficient  depends  primarily  on  wall 
penetration  relative  to  wave  base,  or  relative  to  one-half  of  the  wavelength.  For  instance, 
if  the  bottom  of  the  wall  is  below  the  wave  base,  there  will  be  no  wave  transmission. 

Water  depth  influences  transmission  indirectly  by  influencing  wave  length  in 
intermediate  and  shallow  water.  As  stated  above,  in  deep  water  wave  length  is  dependent 
only  on  period.  The  definition  of  “deep  water”  is  relative  to  water  depth  and  applies 
whenever  the  depth  is  greater  than  or  equal  to  one-half  of  the  wavelength.  As  long  as  the 
wave  base  does  not  touch  the  sea  floor,  the  wave  is  in  deep  water.  As  the  wave  base  begins 
to  come  into  increasingly  greater  contact  with  the  sea  floor,  the  wave  length  is  shortened  due 
to  the  constraint  of  particle  orbits  by  the  sea  floor.  This  shortening  follows  the  equation 


L  = 


tanh  kd 
2tt 


(3) 


where  k  =  2rr/L  and  d  is  the  water  depth.  This  equation  is  applicable  in  intermediate  water, 
which  extends  from  depths  of  L/2  to  1720.  Finally,  as  the  wave  approaches  the  shoreline  and 
passes  a  depth  of  L/20,  it  enters  shallow  water  and  its  length  is  given  by 
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L  =  Tfid  (4) 

Wall  penetration  (w)  is  a  crucial  factor  because  it  can  be  adjusted  by  the  designer  to  control 
the  amount  of  transmission.  Past  wave  records  can  be  analyzed  and  the  data  can  then  be 
used  to  predict  the  wave  periods  that  the  vertical  wall  breakwater  will  be  subject  to.  The 
maximum  expected  wave  periods  and  therefore  wave  lengths  will  have  to  be  considered  in 
designing  a  vertical  wave  barrier.  If  a  designer  wishes  to  eliminate  all  wave  action  behind 
the  barrier,  he  must  only  extend  the  barrier  past  the  wave  base.  Normally,  however,  a  certain 
amount  of  transmitted  energy  is  permitted  in  the  design  in  the  interest  of  monetary  savings, 
and  this  makes  the  determination  of  the  wall  penetration  a  critical,  site-specific,  and 
potentially  risky  decision.  As  will  be  shown,  reducing  wave  barrier  penetration  leads  to  a 
reduction  in  the  wave  force  on  the  barrier,  although  this  is  at  the  expense  of  increased  wave 
transmission. 

Wave  reflection  also  plays  an  important  role  in  determining  wave  forces.  When  a 
wave  reflects  from  a  vertical  wall,  it  no  longer  appears  to  propagate  or  move  in  the  direction 
of  travel.  Rather,  the  wave  seems  to  be  fixed  in  place  with  the  water  surface  moving  only 
up  and  down.  This  condition  is  referred  to  as  a  “standing  wave.”  By  summing  the  equations 
for  water  particle  motion  and  crest  elevation  of  the  incident  and  reflected  waves,  the  motion 
of  the  water  particles  and  wave  crest  elevation  in  a  standing  wave  can  be  determined.  It 
must  be  remembered  that  the  reflected  wave  is  moving  in  the  opposite  direction  of  the 
incident  wave  and  thus  constructive  and  destructive  interference  will  occur  when  these 
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waves  are  superimposed.  Thus,  in  a  standing  wave,  at  the  wall  location,  or  antinode,  all 
water  particle  motion  is  vertical  and  the  surface  elevation  is  equal  to  the  sum  of  incident  and 
reflected  wave  heights.  The  maximum  crest  amplitude  at  the  wall  is  equal  to  (Hj  +  FQ/2  and 
is  repeated  every  one-half  wave  length  away  from  the  wall,  as  shown  in  Figure  6.  At  the 
antinode,  or  trough,  particle  motion  is  horizontal  and  the  surface  elevation  is  equal  to  the 
difference  of  incident  and  reflected  wave  heights,  or  (Hj  -  H,.)/2. 

The  coefficient  of  reflection  is  a  measure  of  the  amount  of  reflection  and  is  given  by 
the  equation 


(5) 


In  a  perfectly  reflecting  wave,  this  coefficient  will  equal  unity,  and  reflected  height  will  equal 
incident  height.  For  vertical  wave  barriers,  reflection  will  be  incomplete  and  Kj  will  be  less 
than  one  since  some  energy  is  transmitted  under  the  wall  while  other  energy  may,  in  fact,  be 
dissipated. 

Incident  wave  height  is  important  for  two  different  reasons.  First,  it  has  a  great 
influence  on  determining  the  amount  of  energy  in  a  wave,  given  by  the  relation 


E  =  Ip gH'L 

o 


(6) 


where  p  denotes  the  specific  gravity  of  water  and  H  is  wave  height.  Second,  if  a  wave  is 
sufficiently  high,  a  portion  of  it  may  pass  over  the  top  of  the  breakwater.  This  is  called 
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overtopping  and  may  result  in  the  transmission  of  a  large  amount  of  energy.  Since  most 
wave  barriers  are  designed  to  be  sufficiently  high  to  prevent  overtopping,  however, 
conservation  of  wave  energy  over  a  flat  bottom  requires  that  all  incident  wave  energy  must 
pass  underneath,  be  reflected,  or  be  dissipated  by  the  barrier.  This  is  shown  in  the  equation 

Ei  =  Et+  Er  (7) 

where  E  is  the  energy  density  and  the  subscripts  I,  r,  and  t  refer  to  the  incident,  reflected,  and 
transmitted  waves  respectively.  When  equation  (6)  is  substituted  into  equation  (7), 
conservation  of  energy  for  partially  transmitted,  partially  reflected  waves  gives 

H?  =  H2r  (8) 

and  thus 

1  (9) 

As  a  wave  impacts  a  full-depth  wall,  all  wave  energy  must  be  reflected.  The 
amplitude  of  the  reflected  wave  will  add  to  the  amplitude  of  the  incident  wave  when  a  wave 
crest  is  at  the  wall.  This  increased  wave  height  will  increase  the  dynamic  pressure  on  the 
wall.  The  force  on  a  wall  due  to  a  perfect  standing  wave  is  then  given  by  the  equation 
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/  =  2&1  +  tanh  kd  +  (10) 

2  k  2 

where  f  is  force  per  unit  length,  the  first  term  computes  the  hydrostatic  pressure  on  the  wall 
and  the  second  term  computes  the  force  due  to  dynamic  pressures  of  the  wave  acting  below 
the  still  water  level.  The  third  term  computes  the  force  due  to  dynamic  pressure  in  the  crest 
of  the  wave  above  the  still  water  level.  If  this  equation  is  applied  to  a  wall  with  water  on 
both  sides,  the  first  term  will  be  balanced  by  the  hydrostatic  pressure  from  water  on  the  other 
side  of  the  wall.  The  large  force  from  the  first  term  that  results  from  dynamic  pressure  below 
the  still  water  level  is  called  a  “first  order  force.”  The  smaller  dynamic  pressure  on  the  wall 
due  to  the  crest  of  the  wave  is  the  source  of  the  “second  order  force.” 
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LITERATURE  REVIEW 


WAVE  TRANSMISSION 

As  previously  mentioned,  studies  on  vertical  wave  barriers  have  been  limited. 
Professor  Robert  Wiegel  of  the  University  of  California  at  Berkeley  is  responsible  for  most 
of  the  theory  that  has  been  developed  for  application  to  vertical  wall  breakwaters.  Wiegel 
(1960)  used  two  different  theories  to  predict  the  ratio  of  transmitted  to  incident  wave 
heights,  or  reflection  coefficient.  A  deep  water  theory  was  developed  by  Ursell  in  1947  and 
states  that 


Ht  _  Kfinw/L) 

H‘  \Itz2I?(2tiw/L)  +K?(2nw/L) 


(11) 


where  w  is  the  wall  penetration  and  K,  (27iw/L)  and  ,1  (2ttw/L)  are  modified  Bessel 
functions.  Wiegel  then  developed  a  theory  for  wave  transmission  valid  in  all  depths  through 
the  use  and  manipulation  of  power  transmission  theory.  Starting  with  the  equations  for 
dynamic  pressure  and  water  particle  velocity,  Wiegel  determined  the  power  transmitted  past 
a  vertical  plane  based  on  linear  wave  theory.  After  some  mathematical  manipulation  of 
terms,  he  determined  that 


K,  = 


H, 

H, 


2 k(d-w)  +  sinh  (2kd-2kw) 
2 kd  +  sinh  2  kd 


(12) 
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In  laboratory  experiments  Wiegel  determined  that  the  power  transmission  theory 
predicted  laboratory  results  more  accurately  than  the  deep  water  theory  for  small  values  of 
barrier  penetration.  This  was  verified  by  Mattson  and  Cederwall  (1976)  as  shown  in  Figure 
7.  Wiegel  found  the  opposite  to  be  true  for  larger  barrier  penetrations.  He  also  observed 
a  trend  of  decreasing  transmission  with  increasing  wave  steepness.  His  final  conclusion  was 
that  wave  barriers  should  be  used  with  caution  because  they  don’t  significantly  prevent  wave 
transmission  unless  they  have  a  large  penetration. 

Mattson  and  Cederwall  (1976)  performed  a  series  of  experiments  to  measure  the 
coefficient  of  transmission  under  various  conditions  and  then  compared  these  measurements 
to  Wiegel  theory.  Their  findings  verified  Wiegel’s  conclusion  that  predictions  for 
transmission  varied  with  barrier  penetration.  This  is  shown  in  Figure  8.  They  also  concluded 
that  Wiegel’s  theory  is  relatively  inaccurate  in  extreme  cases.  Mattson  and  Cederwall  also 
concluded  that  transmission  is  only  a  function  of  wall  penetration  and  wavelength.  They  did 
this  by  plotting  their  experimental  coefficients  of  transmission  versus  various  non- 
dimensional  parameters.  Penetration  divided  by  wavelength  (w/L)  is  the  only  parameter  that 
resulted  in  reasonable  data  correlation.  Through  curve  filtering  they  developed  a  new 
equation  for  the  coefficient  of  transmission,  based  solely  on  incident  wavelength  and  wall 
penetration: 


*,  =  [ 


cosh  O.SAw 


] 


cosh  kw 


(13) 


Figure  7.  Wiegel  theory  compared  to  measured  values  of  K, ,  Mattson  and  Cederwall 
(1976),  pp.  7. 


Figure  8.  Graph  of  the  dependence  of  Kt  on  penetration  relative  to  wave  length  (w/L), 
Mattson  and  Cederwall  (1976),  pp.  7. 
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Another  interesting  finding  was  that  irregular  waves  appeared  to  have  a  greater  coefficient 
of  transmission  than  regular  waves.  Mattson  and  Cederwall  concluded  that  deviations  from 
Wiegel  theory  was  most  likely  the  result  of  energy  loss  due  from  vortex  formation  as  water 
flowed  underneath  the  breakwater.  Also,  wave  reflection  from  the  wall  generates  non-linear 
conditions  in  front  of  the  wall  which  may  invalidate  wave  linearity  assumptions  made  by 
Wiegel. 

Losada  et  al  (1992)  examined  the  applicability  of  linear  theory  to  wave  interactions 
with  vertical  wall  breakwaters.  An  eigenfunction  expansion  method  was  applied  to  obtain 
theoretical  solutions  for  transmission  and  reflection  coefficients.  This  paper  utilized  two 
assumptions  that  were  adopted  in  the  theoretical  development  of  this  Trident  project 
concerning  velocity  and  pressure  at  the  breakwater:  Both  pressure  and  velocity  were 
assumed  to  be  constant  across  the  gap  underneath  the  barrier.  The  numerical  analysis 
demonstrated  in  this  paper  created  a  numerical  theory  which  predicted  much  larger 
coefficients  of  transmission  than  those  from  Wiegel.  Losada  et  al  also  showed  that  forces 
on  the  barrier  were  a  function  of  barrier  penetration  as  well  as  water  depth  and  could  be 
computed  by  vertically  integrating  the  pressure  difference  across  the  wall  using  the  equation 

o 

Az)  =  I  (P2~Px)dz  (14) 

-W 

where  f  is  the  force  on  the  wall,  p2  is  the  pressure  on  the  front  of  the  wall,  pj  is  the  pressure 
on  the  back  of  the  wall,  and  w  is  the  penetration  of  the  wall. 

Several  papers  have  been  written  by  engineers  at  the  company  of  Peratrovich, 
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Nottingham  and  Drage,  Inc.  (PN&D)  regarding  vertical  wall  breakwaters.  They  review  the 
history  and  development  of  these  breakwaters,  describe  tests  conducted  on  breakwater 
models  for  feasibility  studies,  and  conclude  that  vertical  wall  breakwaters  are  very  effective 
at  decreasing  wave  transmission  and  perform  as  designed.  One  report  notes  that  in  the  case 
of  a  vertical  wall  breakwater  in  Garibaldi,  Oregon,  the  structure  stood  undamaged  from 
storm  conditions  twice  as  severe  as  design  conditions  (PN&D,  Issue  #30,  undated).  It  is  not 
clear,  however,  whether  these  observations  were  made  regarding  incident  wave  heights  or, 
more  likely,  the  sum  of  incident  and  reflected  heights. 

Gilman  and  Shaver  (1992)  reported  on  tests  performed  on  a  model  of  a  PN&D  wave 
barrier.  The  tests  measured  wave  transmission  and  loadings  in  hopes  of  optimizing  designs 
for  proposed  wave  barrier  installations  in  Seward,  Alaska  and  Seattle,  Washington.  The 
model  was  constructed  of  epoxy-coated  sheets  of  plywood  connected  to  an  aluminum  and 
steel  frame  which  allowed  elevation  and  slope  changes  of  the  barrier.  Geometric  scaling  was 
at  a  1:12  ratio  with  velocities  and  times  scaled  according  to  Froude  modeling  laws.  The 
model  was  subjected  to  regular  waves  with  full-scale  periods  between  three  and  six  seconds 
and  irregular  wave  trains  based  on  JONS  WAP  frequency  spectra  in  a  full-scale  water  depth 
of  18.3  meters.  During  testing,  the  barrier  limited  transmission  to  25%  when  penetration-to- 
depth  ratios  were  greater  than  approximately  0.3.  Gilman  and  Shaver  (1992)  noted  that 
most  deep  water  waves  were  completely  reflected  when  penetration  was  greater  than  6 
meters  and  that  much  of  the  lower-frequency  wave  energy  was  reflected  in  that  case  as  well. 
They  concluded  that  the  tested  wave  barrier  could  withstand  wave  heights  of  three  meters 
in  real  life  and  would  therefore  be  very  useful  in  areas  not  completely  subjected  to  open- 
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ocean  wave  assault. 


EQRCES  ON  VERTICAL  WAVE  BARRIERS 

The  Army  Corps  of  Engineers’s  Shore  Protection  Manual  (1984)  and  Naval 
Facilities  Engineering  Command  Design  Manual  26.2  (1982)  contain  an  extensive  analysis 
of  wave  forces  on  vertical  walls.  Although  not  specific  to  wave  barriers,  the  theory  is  very 
similar  and  can  be  applied  to  vertical  wave  barriers  with  a  few  simple  modifications.  The 
Shore  Protection  Manual  (SPM)  relies  on  the  aforementioned  Miche-Rundgren  method  of 
determining  non-breaking  wave  forces  (Figure  9).  In  this  method,  the  height  of  water  at  the 
wall  is  determined  to  be  equal  to  the  sum  of  incident  and  reflected  wave  heights  and  can  be 
assumed  to  be  twice  the  incident  wave  height  with  perfect  reflection  (K,  =  1).  The  water 
height  will  oscillate  about  a  point  (h0)  above  the  still  water  level.  This  elevation  is  located 
midway  between  the  wave  crest  and  wave  trough  and  is  greater  than  zero  for  nonlinear 
waves  because  such  waves  are  asymmetric  with  crest  amplitudes  that  are  greater  than  trough 
amplitudes.  Wave  forces  on  the  wall  are  determined  to  be  the  sum  of  hydrostatic  (ph  =  pgd) 
and  dynamic  pressures  (Figure  10).  The  dynamic  pressure  at  the  sea  floor  is  given  by: 

_  (l+Kr)  pgn, 

Pd  1  2  cosh  (2  nd/L) 


(15) 


Crest 
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Figure  9.  Sketch  of  terms  in  Miche-Rundgren  method  of  predicting  forces  on  walls.  U.S. 
Army  Corps  of  Engineers  (1984),  pp.  7-162. 


Figure  10.  Hydrostatic  and  dynamic  pressures  on  a  wall.  U.S.  Army  Corps  of  Engineers 
(1984),  pp.  7-163. 
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When  a  wave  crest  is  at  the  wall,  the  force  on  the  wall  per  unit  width  is  then  approximated 
by 


/=  |(P gd+Pd){d  +  hQ  +H) 


(16) 


This  force  is  then  presented  in  a  dimensionless  form  given  by 


pgd2  2 


(17) 


Similar  expressions  can  be  derived  for  the  force  when  a  wave  trough  is  at  the  wall.  Example 
graphs  from  the  SPM  for  determining  these  forces  are  shown  in  Figures  1 1  and  12.  These 
graphs  require  several  secondaiy  calculations  and  can  cause  error  in  force  prediction  because 
of  the  interpolation  that  is  necessary  to  read  data. 

The  SPM  then  develops  theory  and  contains  an  example  for  the  force  on  a  wall  built 
on  a  rubble  foundation.  As  shown  in  Figure  13,  the  distribution  of  force  in  this  case  is 
extremely  similar  to  that  on  a  vertical  wave  barrier.  The  force  in  this  case  is  shown  to  be 


f"  =  (l-i»/ 


(18) 


where  rf  is  simply  a  force  reduction  factor 


rr 


d-w  _  d-w 
d  (  ~d~ 


(19) 


Figure  11.  Required  figure  for  force  prediction  using  Miche-Rundgren  method.  U.S.  Army 


Corps  of  Engineers  (1984),  pp.  7-164. 


Figure  12.  Required  figure  for  force  prediction  using  Miche-Rundgren  method.  U.S.  Army 
Corps  of  Engineers  (1984),  pp.  7-165. 
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Figure  13.  Distribution  of  force  on  a  wall  atop  rubble  foundation,  a  case  similar  to  a  wave 
barrier.  U.S.  Army  Corps  of  Engineers  (1984),  pp.  7-178. 
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proportional  to  the  height  of  the  foundation,  d  -  w,  divided  by  the  total  height  of  the  wall, 
d. 

The  NAVFAC  Engineering  Design  Manual  (1982)  contains  an  example  problem 
specific  to  vertical  wave  barriers.  The  equations  for  computing  forces  on  the  barrier  are  the 
same  as  those  used  by  the  SPM  in  the  rubble  foundation  example.  This  example 
problem  and  method  were  used  as  examples  to  which  to  compare  project-generated 
experimental  and  theoretical  data.  One  critical  assumption  employed  by  both  the  NAVFAC 
and  SPM  methods  is  that  a  wave  crest  is  on  the  incident  side  of  the  barrier  while  a  wave 
trough  is  on  the  leeward  side  of  the  barrier.  As  noted,  this  produces  the  maximum  possible 
pressure  difference  across  the  wall  and  thus  the  maximum  force  on  the  wall.  Such  an 
occurrence,  however,  seems  unrealistic  and  appears  to  predict  unrealistically  large  forces  on 
a  wall.  As  will  be  shown  in  later  sections,  these  forces  are  indeed  unrealistically  large, 
typically  twice  as  large  as  those  measured  in  laboratory  experiments. 


RESEARCH  OBJECTIVES 

Based  on  the  review  of  existing  literature,  it  was  determined  that  there  was  a  need 
for  additional  research  in  the  areas  of  wave  transmission  and  wave  forces.  Specifically,  this 
additional  research  in  the  area  of  wave  transmission  included  the  development  of  a  modified 
theory  of  wave  transmission,  a  comparison  of  wave  transmission  data  from  the  laboratory 
to  available  theories,  and  an  evaluation  of  the  modified  theory  for  wave  transmission  when 
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compared  to  actual  test  results.  Additional  research  in  the  area  of  wave  forces  was  divided 
into  the  same  groups  as  wave  transmission:  Development  of  theory,  comparison  of  data  to 
conventional  theory,  and  evaluation  of  developed  theory. 

As  a  result,  two  main  thrusts  were  developed.  The  first  was  to  develop  a  modified 
theory  for  wave  transmission  and  wave  forces.  The  second  was  to  then  test  these  theories 
in  the  laboratory  by  comparing  experimental  data  to  predictions  from  the  newly  developed 


theories. 
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THEORY 


POWER  TRANSMISSION 

The  modified  theory  for  wave  transmission  and  wave  force  on  a  vertical  wave  barrier 
is  an  extension  of  existing  theory  as  developed  by  Wiegel  (1960).  Wiegel  proposed  that 
wave  transmission  could  be  predicted  from  the  wave  power  passing  under  a  wave  barrier. 
To  make  this  prediction,  Wiegel  assumed  that  there  was  no  energy  loss  in  the  system;  all 
energy  was  either  transmitted  beneath  the  barrier  or  reflected.  According  to  Wiegel,  if  a 
vertical  power  distribution  is  considered  such  as  in  Figure  14,  that  portion  of  the  wave  power 
that  is  between  the  bottom  of  the  barrier  and  the  sea  floor  will  be  transmitted.  All  other 
incident  power  will  be  reflected. 

Wave  power  per  unit  area  is  the  product  of  the  dynamic  pressure  (force  per  unit 
area)  and  horizontal  particle  velocity  at  any  point  in  a  wave.  Thus  the  total  power  present 
in  a  complete  wave  passing  through  an  imaginary  vertical  section  in  two-dimensions  is  given 
by  the  time-averaged  product  of  wave-induced  dynamic  pressure  and  particle  velocity  at  this 
section,  integrated  from  the  tank  floor  to  the  still-water  level.  The  initial  power  without  any 
wave  barrier  present  is 


P  = 


I  j  \pudzdt 

1  0  -d 


(20) 
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where  P  =  wave  power  p  =  pressure  from  linear  wave  theory,  and  u  =  the  horizontal  particle 
velocity  as  given  in  the  equations 


P  = 


pgH  cosh  k(z+d) 
2  cosh  kd 


cos  ( kx-ot ) 


(21) 


_  oi¥  cosh  k(d+z) 
2  sinh  kd 


cos  (kx-ot) 


(22) 


where  o  is  the  angular  frequency  of  the  wave  and  is  equal  to  2rc/T.  When  integrated,  (20) 
yields 


P  -Ipgff’C, 


(22) 


where  l/8pgH2  is  the  wave  energy  density  and  where  Cg  is  the  wave  group  velocity  given  by 


C  „  2kd  . 
—  (1+ - ) 

2  sinh  2  kd 


(23) 


In  equation  (23),  C  is  the  wave  celerity  and  is  equal  to  L/T.  Equation  (22)  is  important 
because  it  states  that  wave  power  in  a  closed  system  is  constant,  independent  of  time  or 
location,  when  measured  over  one  complete  wave  cycle. 
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WIEGEL  THEORY 


As  stated  earlier,  according  to  Wiegel  power  will  only  be  transmitted  across  the 
section  of  the  tank  not  blocked  by  the  wave  barrier.  In  other  words,  power  will  only  be 
transmitted  across  the  gap  beneath  the  bottom  of  the  wall.  Therefore,  the  product  of 
incident  wave  pressure  and  velocity  (equation  20)  need  only  be  integrated  between  the  sea 
floor  (-d)  and  the  bottom  of  the  wall  (-w).  This  generates  the  following  equation  for  power 
transmission  past  a  wall: 


P,  =  J  J  P,  u,  d2  dt 

1  0  -d 

Which,  after  integration,  gives 


(24) 


—osH2C  r  Ikid-w) +sinh(2fc/-2£w) , 
8  '  8  2fo/+sinh  2  kd 


(25) 


To  find  the  transmission  coefficient,  Wiegel  reasoned  that  the  power  downstream  (behind 
the  wall)  must  equal  the  power  transmitted  beneath  the  wall.  Considering  a  vertical  section 
downstream. 


p,  =  7  J  lPtu,dzdt 

1  0  -d 


(26) 


and  the  time-averaged  transmitted  power  is  given  by 
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(27) 


Wiegel  determined  the  coefficient  of  transmission  (as  stated  earlier  in  Equation  1)  by 
equating  (25)  and  (27)  such  that 


Ht 

K =-i  = 
'  77. 


1 


2  k(d-w)  +  sinh(2kd-2kw) 
Ikd  +  sinh  2 kd 


(28) 


MODIFIED  TRANSMISSION  THEORY 

When  deriving  the  power  transmission  underneath  the  wave  barrier,  Wiegel  did  not 
consider  the  effects  of  reflection  upon  wave  transmission.  Consider  a  lone  wave  approaching 
a  vertical  wave  barrier.  As  the  crest  of  the  wave  begins  to  impact  the  wall,  the  water  level 
at  the  wall  is  rising.  As  this  is  happening,  part  of  the  wave’s  energy  is  transmitted  beneath 
the  wall  and  continues  progressing  down  the  tank.  The  part  of  the  wave’s  energy  that  is 
reflected  begins  to  travel  back  towards  the  wave  maker.  This  reflected  wave  will  then 
interfere,  either  constructively  or  destructively,  with  the  incident  wave  energy  that  is  still 
approaching  the  wall. 

When  the  highest  point  of  the  wave  crest  is  at  the  wall,  the  water  level  is  at  its 
maximum  value  because  the  incident  and  reflected  crests  constructively  interfere.  With 
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perfect  reflection  the  water  level  at  the  wall  will  be  equal  to  twice  the  incident  wave  height. 
Dynamic  pressure  will  also  be  at  its  maximum  and  equal  to  the  sum  of  incident  and  reflected 
pressure.  Water  particle  velocities,  however,  are  not  at  their  maximum  value  because  the 
incident  and  reflective  wave  velocities  destructively  interfere  with  each  other.  The  horizontal 
velocities  of  the  incident  wave  particles  are  towards  the  wall  while  the  particles  of  the 
reflected  wave  are  moving  away  from  the  wall.  Thus  when  a  wave  crest  is  at  the  wall, 
pressures  are  additive  (total  pressure  equals  incident  plus  reflected)  due  to  constructive 
interference  and  as  a  result  of  destructive  interference  velocities  are  subtractive  (net  velocity 
equals  incident  minus  reflected). 

The  development  of  a  modified  theory  for  wave  transmission  begins  with  a 
modification  of  the  initial  power  transmission  equation  to  account  for  reflection,  as  was 
already  discussed.  In  the  case  of  a  vertical  wave  barrier,  where  reflection  is  not  complete, 
pressures  are  additive  but  are  not  equal  to  twice  the  incident  pressures  and  velocities  are 
subtractive  but  do  not  completely  cancel  each  other.  The  transmitted  power  can  then  be 
found  by  integrating  net  pressures  across  the  wall  and  net  particle  velocities  across  the  gap 
beneath  the  barrier  such  that 


~W 


I  (Pi+PXUrUr)  &  dt 

-d 


(29) 


Linear  theory  is  used  to  provide  approximations  for  pressures  and  velocities.  After 
substitution  and  integration,  the  final  modified  equation  for  power  transmitted  under  a  wall 
is 
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(30) 


Some  interesting  comparisons  can  be  made  between  this  equation  and  the  wave 
power  equation  without  the  effects  of  reflection  (25).  From  (30)  it  is  evident  that  any 
reflection  will  decrease  the  modified  prediction  for  net  power  transmission  from  that  of 
Wiegel’s  (equation  25).  When  there  is  complete  reflection  (K,  =  1)  there  is  no  power 
transmission.  When  there  is  no  wall  (penetration  w  =  0  and  hence  no  reflection),  the 
reflection  term  as  well  as  the  last  term  (in  brackets)  go  to  unity  and  the  original  equation  for 
net  time-averaged  power  transmission  without  a  wall  (22)  remains.  Under  any  conditions 
where  there  is  a  wall,  however,  the  last  term  associated  with  wall  penetration  is  less  than 
one,  the  reflection  coefficient  is  greater  than  zero,  and  the  power  transmitted  is  less  than  that 
predicted  by  Wiegel. 

When  the  modified  wave  power  equation  (30)  is  equated  to  transmitted  power 
downstream  from  the  wave  barrier  in  equation  (27)  ,  the  coefficient  of  transmission  can  be 
determined  by  simplification  to  be 


K]  =  (!-*’)[ 


2£(<7-w)+sinh  2k(d-w) 
2kd+sm\\  2  kd 


(31) 


Further  simplification  is  possible  if  conservation  of  momentum  (mass  flow)  beneath  the  wall 
is  assumed.  It  is  recognized  that  fluid  velocities  under  the  barrier  are  continuous  so  that 
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w,  =u-ur  (32) 

After  substitution  of  linear  wave  theory  terms  for  particle  velocity  and  simplification,  it  is 
found  that 

K,  =  1  -Kr  (33) 

Finally,  combining  equations  (31)  and  (33)  gives  the  desired  solution  for  the  transmission 
coefficient 


K  _  £  ri+  2fo/  +  sinh  2kd  -..i 
'  1  2k(d-w)  +  sinh  2k(d-wy 


(34) 


The  assumption  of  conservation  of  momentum  instead  of  conservation  of  energy 
allows  for  a  loss  of  energy  in  the  system.  Assuming  conservation  of  flow  requires  a 
modification  of  equation  (7)  to  account  for  the  lost  energy  such  that 


E{  -  Er+  E  +  Ewss 


(35) 


If  energy  loss  is  considered  relative  to  the  incident  wave  energy  and  simplifications  of  the 
resulting  equation  are  made,  a  relation  for  energy  loss  is  obtained  such  that 


■'LOSS 


(36) 
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This  equation  shows  that  the  maximum  energy  loss  is  one-half  of  the  incident  energy  and 
will  occur  when  the  coefficient  of  transmission  equals  the  coefficient  of  reflection,  which 
will  only  happen  when  K,  =  K,  =  0.5. 


MODIFIED  FORCE  THEORY 


The  force  on  a  vertical  wave  barrier  can  be  determined  by  the  integration  of  the 
pressure  difference  between  the  front  and  back  of  the  wall.  The  front  of  the  wall  is 
subject  to  the  sum  of  incident  and  reflected  pressures  while  the  rear  of  the  wall  is  only 
subjected  to  transmitted  wave  pressure.  Thus,  the  net  pressure  on  the  wall  at  any  instant 
in  time  equals  the  pressure  at  the  front  of  the  wall  minus  the  pressure  at  the  back  of  the 
wall  or  Pj  +  Pr  -  pt.  When  these  pressures  are  to  be  integrated  over  the  height  of  the  wall, 
thereby  yielding  force  per  unit  width  on  the  wall,  the  equation  becomes 

f=Hpi+pr-p,)dz  <37> 

-w 

When  this  equation  is  integrated,  accounting  for  a  90  degree  phase  shift  between  the 
incident  and  transmitted  waves,  it  is  found  that 


P gHt 

2k 


(cos  at  (1  +.£.)- sin  at  (Kt )) 


[sinh  M-sinh  k(d-w 
cosh  kd 


(38) 


The 
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modified  theory  was  initially  derived  with  the  assumption  that  water  levels  on  both  sides 
of  the  wall  would  be  in  phase.  This  contradicts  the  180  degree  phase  shift  assumed  in 
NAVFAC  calculations,  but  was  believed  to  be  more  reasonable.  After  the  first  test 
results  were  analyzed,  though,  it  was  discovered  that  there  is  a  90  degree  phase  shift  in 
water  levels  between  the  front  and  back  of  the  wave  barrier,  as  shown  in  Figure  15. 
Therefore  when  a  wave  crest  is  at  the  wall,  the  water  level  behind  the  wall  will  be  at 
approximately  the  still  water  level.  One-half  second  later  (for  a  0.5  Hertz  wave  such  as  in 
Figure  15)  the  crest  will  be  at  the  back  side  of  the  wall.  When  the  maximum  value  of  the 
sine-cosine  term  is  used  to  predict  forces,  the  previous  equation  can  be  further  simplified 
using  this  assumption  and  through  application  of  equation  (33)  to  give 

F. „ KH,  (39) 

The  above  equation  accounts  for  the  dynamic  pressures  due  to  waves,  or,  when 
integrated,  the  first  order  forces.  First  order  forces  are  the  dominant  factor  in  predicting 
the  total  force  on  the  wall.  In  the  free  surface  region  between  the  still  water  level  and  the 
wave  crest,  dynamic  pressures  are  assumed  to  be  like  hydrostatic  pressures:  Increasing 
from  pd  =  0  at  the  crest  to  pd  =  pgh  at  the  still  water  level.  The  difference  in  these 
pressures  due  to  the  crest  can  also  be  integrated  across  the  wave  barrier  to  obtain  another 
expression  for  force.  This  second  expression  for  force  is  called  the  second  order  force. 
This  force  is  not  as  large  as  the  first  order  force,  but  can  account  for  a  significant  amount 
of  the  total  force  when  conditions  with  large  wave  heights  are  examined.  The  equation 


Figure  15.  Graph  of  watei 
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for  this  second  order  force  is 


8 


(40) 


where  fc  is  the  adjustment  for  the  crest  at  the  front  of  the  wall  and  H*.  is  the  water  height 
at  the  wall.  If  this  adjustment  is  simplified  using  equation  (33)  and  added  to  equation 
(39)  the  final  modified  equation  for  the  maximum  dynamic  force  on  a  vertical  wave 
barrier  is 


/ 


_  fagH,  J2_2k  +£2  [sinh  fc/-sinh  kjd-w)]  +  PgH,  (2  K) 
h  '  11  cosh  kd  fi 


(41) 


This  is  a  very  important  equation  because  it  specifies  which  wave  and  environmental 
characteristics  determine  the  force  on  a  vertical  wave  barrier.  From  this  equation  it  is 
apparent  that  force  depends  on  wave  amplitude,  wave  length,  water  depth,  and  wall 
penetration.  Since  wave  length  is  in  fact  determined  by  wave  period  and  water  depth,  a 
functional  dependence  on  wave  period  is  also  indicated.  One  key  feature  of  equation  (41) 
is  that  force  depends  on  K,  as  calculated  from  equation  (34).  This  presents  a  challenging 
cost-benefit  analysis  problem  because  while  the  designer  of  a  wave  barrier  may  seek  to 
minimize  wave  transmission  past  a  wave  barrier  into  a  harbor,  this  can  only  be  achieved 
at  the  expense  of  larger  forces  on  the  barrier.  Larger  forces  then  require  the  construction 
of  a  larger,  stronger,  more  expensive  breakwater  that  can  withstand  these  increased 


forces. 
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EXPERIMENTAL  SETUP 


GENERAL  SETUP 


Testing  of  the  laboratory  model  was  conducted  in  the  United  States  Naval 
Academy  Hydromechanics  Laboratory’s  120-foot  wave/towing  tank.  Overall  tank 
dimensions  are  120  feet  in  length,  8  feet  in  width,  and  5  feet  in  depth  (Figure  16).  The 
tank  is  fabricated  from  concrete  except  for  a  25  foot  segment  from  40  to  65  feet 
(measured  from  the  wave  maker)  where  the  starboard  tank  wall  is  constructed  of  glass 
viewing  windows.  A  dual-flap  type  wave  maker  is  used  to  generate  regular  and  random 
waves  with  frequencies  between  0.2  and  1.3  Hz  and  heights  up  to  approximately  8 
inches.  These  waves  are  then  dissipated  at  the  other  end  of  the  tank  by  a  sloped  beach 
consisting  of  lattice  bundles  of  fiberglass  rods  and  stainless  steel  connectors. 

The  model  vertical  wave  barrier  was  placed  a  little  more  than  60  feet  from  the 
wave  maker  and  was  thus  slightly  closer  to  the  beach  than  the  wave  maker  because  the 
presence  of  the  glass  wall  section  precluded  the  mounting  of  the  wall  at  the  exact 
midpoint  of  the  tank.  This  allowed  valid  data  to  be  taken  for  15  to  30  seconds, 
depending  on  the  generated  wave,  before  waves  reflected  from  the  beach  end  of  the  tank 
returned  to  the  wall,  thus  contaminating  the  data.  For  tests  with  regular  waves,  incident 
waves  were  measured  with  a  gage  15  feet  from  the  wave  maker.  Transmitted  waves 
were  measured  with  a  gage  75  feet  from  the  wave  maker,  or  15  feet  behind  the  barrier. 
Water  levels  at  the  barrier  were  measured  using  one  wave  gage  one  inch  in  front  of  the 


IPO'  Vave/Towing  Tank  --  Plan  View 
Regular  Wave  Testing  Conf  iguration 


Figure  16.  Plan  view  of  wave/towing  tank  for  regular  wave  setup. 


barrier  and  another  gage  two  inches  behind  the  barrier.  The  incident  and  transmitted 
wave  gages  were  resistance  wire  probes  while  those  against  the  wall  were  capacitance 
probes.  Additionally,  a  sonic  wave  probe  was  mounted  on  a  carriage  that  runs  above  the 
tank  for  most  of  its  length.  The  sonic  gage  (carriage)  began  its  run  above  the  model  and 
moved  in  the  direction  of  the  wave  maker  in  order  to  record  a  water  surface  profile 
during  the  experiment  that  could  be  analyzed  to  determine  incident  and  reflected  wave 
heights. 


SCALING  LAWS 


The  entire  experiment  was  modeled  on  a  1:10  scale  in  order  to  produce  results 
that  could  be  transferred  to  a  real-life  50  foot  water  depth  according  to  geometric  scaling. 
In  this  case,  the  modeling  of  a  50-foot  water  depth  by  a  5-foot  deep  tank  resulted  in  a 
length  ratio  of  10  such  that 


(42) 


where  Lr  is  the  length  ratio,  Lp  is  the  prototype  length,  and  Lm  is  the  model  length. 
Froude  scaling  is  used  in  this  problem  as  Froude  scaling  is  utilized  whenever  surface 
gravity  wave  motions  are  to  be  simulated.  The  equation  for  the  Froude  number  is 
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where  U  is  a  characteristic  velocity,  L  is  a  characteristic  length,  and  g  is  the  acceleration 
of  gravity.  A  model  and  prototype,  or  real-life  structure,  will  have  the  same  Froude 
number  when  velocity  and  length,  the  two  variables  of  a  Froude  number,  are 
proportionally  related  to  give  the  same  value  of  Fr .  With  equal  Froude  numbers 


In  this  experiment,  the  use  of  a  length  ratio  of  10  resulted  in  a  time  ratio  of  (10)y\ 
Therefore  a  2.0  second  period  wave  from  this  experiment  corresponds  to  a  wave  with  a 
period  of  2(10),/‘  =  6.3  seconds  in  real  life.  Laboratory  test  conditions  ranged  from  0.9 
second  period  to  2.0  second  period  waves  because,  when  scaled,  these  wave  periods 
roughly  correspond  to  real-life  periods  common  to  inland  waters  of  the  Pacific 
Northwest.  (PN&D,  Issue  #30,  undated) 
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MODEL  CONSTRUCTION 

The  model  wall  is  constructed  of  G-10  composite  sheets  laminated  onto  a  waffle 
interior  framework.  The  composite  sheets  are  1/4"  thick  in  front  of  the  wall  and  3/8 
thick  on  the  rear  of  the  structure,  with  the  thicker  sheet  placed  on  the  back  to  increase  the 
stiffness  of  the  wall.  The  interior  waffle  pattern  is  composed  of  square  foam  and 
aluminum  bar.  The  foam  runs  vertically  and  horizontally  and  provides  buoyancy  and 
rigidity  while  three  strips  of  aluminum  bar  are  placed  vertically  to  provide  additional 
rigidity  to  the  wall  as  well  as  to  allow  the  mounting  of  the  wall  onto  a  support  structure. 
The  wall  is  free-flooding  as  a  result  of  holes  that  were  drilled  through  all  horizontal 
section  pieces  and  as  a  result  prevents  an  excessive  buoyant  force  from  developing  when 
the  wall  is  placed  in  the  wave  tank  and  induces  in  a  slight  negative  buoyancy  overall  for 
the  structure. 

The  wall  is  mounted  to  a  support  structure  composed  of  cross-pieces  of 
aluminum  angle  and  support  legs  composed  of  aluminum  angle  and  tubing  (Figure  1 7). 
The  aluminum  cross  pieces  are  3"  x  3"  angle  of  1/4"  thickness.  The  bottom  cross-piece  is 
reinforced  with  a  3/8"  aluminum  plate  bolted  to  the  bottom  edge.  This  plate  was  added 
to  increase  the  stiffness  of  the  aluminum  angle  in  order  to  prevent  flexure.  There  are  slots 
cut  in  the  forward  leg  of  angle  to  allow  mounting  of  the  wall  to  the  angle  through  the  use 
of  1/4"  Allen-head  screws  put  into  holes  drilled  and  tapped  in  the  aluminum  stiffeners  of 
the  wall. 

The  support  frame  consists  of  two  cross  pieces  running  across  the  wave  tank, 
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Figure  17.  Plan  view  of  model  wave  barrier. 
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parallel  to  the  tank  floor  (Figure  18).  They  are  attached  to  the  tank  wall  with  one  inch 
diameter  brass  pegs  placed  through  a  hole  in  the  end  of  the  cross-pieces  into  pre-existing 
holes  in  the  tank  walls.  The  cross-pieces  and  wall  are  further  reinforced  against  flexure 
by  the  use  of  structural  aluminum  members  as  “tie  rods”  which  run  from  the  cross-pieces 
to  fasteners  on  the  tank  wall  eight  feet  behind  the  model.  The  bottom  members  are  one 
inch  diameter  tubing  aluminum,  while  the  top  legs  are  composed  of  aluminum  angle  of 
the  same  size  as  the  cross-pieces.  Tubing  was  used  on  the  bottom  of  the  wall  to  decrease 
interference  with  the  flow  behind  the  wall. 

The  model  wave  barrier  consists  of  three  distinct  sections  .  The  two  outside 
sections  are  fixed  to  the  support  structure  by  1/4"  Allen-head  screws  and  5/16"  through- 
bolts.  These  sections  are  spaced  from  the  wall  by  putting  each  screw  through  two 
stacked  1/8"  washers  before  attaching  to  the  wall.  The  middle  section  is  attached  to  the 
support  structure  only  through  two  pivot-force  block  connections  (Figure  18).  These 
connections  allow  the  measurement  of  force  in  the  direction  of  wave  propagation. 

Rows  of  holes  were  drilled  in  the  wall  at  six  inch  intervals  to  allow  movement 
vertically  without  necessitating  a  movement  of  the  support  structure.  It  was  originally 
intended  to  move  the  wall  in  the  tank  by  removing  the  connecting  screws,  repositioning 
the  wall,  and  fastening  it  to  the  supports  once  again.  However,  due  to  extra  hardware 
used  for  spacing  and  reinforcement  and  the  general  awkwardness  of  wall,  it  soon  became 
readily  apparent  that  this  method  would  be  too  time  consuming  to  use.  After  some 
consideration,  the  wall  was  repositioned  by  removing  the  rear  reinforcing  legs  and  the 
brass  pegs  that  attached  the  model  to  the  tank  walls.  The  whole  wall  was  then  elevated 
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or  lowered  and  the  pegs  were  reinserted  in  the  bottom  cross-piece.  The  top  cross-piece 
was  then  removed  and  reattached  at  its  original  level  at  the  top  of  the  tank  wall.  Finally, 
the  reinforcing  legs  were  reattached. 

The  fine  tolerances  required  in  the  construction  of  this  model,  especially  due  to 
the  use  of  two  force  gages,  made  construction  difficult  and  time-consuming.  Use  of  a 
single  force  gage  to  measure  wave  forces  by  suspending  the  center  panel  in  cantilever- 
fashion  would  have  greatly  simplified  design  and  assembly.  This  would  have  precluded 
possible  determination  of  the  center  of  pressure  or  of  wave-induced  moments  on  the  wall, 
however. 


EXPERIMENTAL  CONDITIONS 


The  first  part  of  the  experiment  involved  the  creation  of  regular  waves  of  fixed 
frequency  and  wave  steepness  in  order  to  test  the  modified  theory  for  wave  transmission 
and  resultant  force.  Waves  with  frequencies  of  0.5,  0.7,  0.9,  and  1.1  Hertz  were 
generated  at  wave  steepnesses  of  1/40,  1/20,  and  1/15  where  steepness  is  equal  to  wave 
height  divided  by  wave  length  =  H/L.  The  use  of  fixed  slopes  instead  of  fixed  wave 
heights  as  a  testing  variable  was  the  result  of  the  observed  dependence  of  the 
dimensionless  wave  force,  and  in  particular  non-linear  wave  forces,  on  the  ratio  of  wave 
height  to  wave  length  (H/L)  as  shown  in  the  modified  theory.  All  tests  were  conducted  in 
a  water  depth  of  five  feet  but  with  three  different  wall  penetrations  of  3.0,  2.5,  and  2.0 
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feet.  In  this  first  part,  waves  were  measured  at  15  feet  and  85  feet  from  the  wave  maker 
as  well  as  the  front  and  back  of  the  wall  (Figure  16).  In  addition,  the  forces  on  the  top 
and  bottom  force  gages  were  recorded.  Data  were  recorded  for  a  30  second  duration 
using  a  sampling  frequency  of  24  Hertz. 

The  goal  of  the  second  phase  of  testing  was  to  measure  fluid  velocities  beneath 
and  in  front  of  the  wall  in  a  grid  pattern,  as  shown  in  Figure  19.  This  was  accomplished 
with  the  aid  of  an  acoustic  Doppler  velocimeter  (ADV)  manufactured  by  Sontek,  Inc.  of 
San  Diego,  California  (Figure  20).  This  device  emits  sonic  waves  which  are  reflected  by 
particles  in  the  water.  The  ADV  then  measures  the  Doppler  shift  of  its  emitted  waves 
and  can  determine  particle  velocities.  A  special  rig  was  created  which  allowed  the  probe 
to  be  mounted  parallel  to  the  tank  bottom,  thereby  minimally  interfering  with  fluid  flow. 
Additionally,  this  rig  allowed  for  precise  alteration  of  the  probe’s  position,  which  was 
essential  for  the  second  part  of  this  testing  phase.  A  0.5  Hertz,  1/40  steepness  wave  was 
generated  for  this  part  of  the  experiment. 

The  third  and  final  phase  of  testing  studied  the  effects  of  random  waves  on  the 
wall.  A  computer  program  was  used  to  generate  the  random  waves  in  this  phase.  The 
program  requested  the  entry  of  a  desired  spectrum  shape,  spectrum  width,  modal  (peak) 
period,  significant  height,  and  number  of  frames  and  then  proceeded  to  generate  a  model 
spectrum  by  superimposing  sinusoidal  waves  of  different  frequencies  and  amplitudes  in 
order  to  approximate  a  spectral  shape.  This  spectrum  could  then  be  created  in  the  wave 
tank,  analyzed  by  the  computer,  and  input  into  the  wave  maker  after  being  adjusted  to 
create  a  more  accurately  distributed  spectrum. 
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Figure  20.  Photo  of  an  Acoustic  Doppler  Velocimeter. 
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This  program  was  used  to  create  wave  records  60  seconds  in  length,  of  0.5  and  0.7  Hertz 
modal  frequency,  with  a  JONS  WAP  spectral  form  having  a  spectral  width  of  3.3.  The 
JONSWAP  spectral  equation  is  given  on  page  138  of  Sorenson  (1993).  The  significant 
height,  H, ,  is  a  measure  of  the  wave  height  of  the  1/3  largest  waves  of  a  spectrum.  A 
significant  wave  height  was  chosen  that  produced  a  wave  steepness  H,  /L  =  1/40.  The 
significant  wave  height  was  therefore  5.7  inches  for  a  0.5  Hertz  peak  frequency  and  3.1 
inches  for  the  0.7  Hertz  peak  period.  Five  60-second  time  series  were  created  for  each 
modal  frequency  thus  subjecting  the  model  to  about  100  total  waves  at  0.5  Hertz  and 
approximately  140  waves  at  0.7  Hertz  frequency.  Additionally,  two  wall  penetrations  of 
2.5  and  3  feet  were  tested.  In  total,  four  different  experimental  conditions  were 
examined. 

Force  on  the  wall,  incident  and  transmitted  wave  heights,  and  water  level  at  the 
front  and  rear  of  the  wall  were  the  variables  measured  in  this  section  of  the  experiment. 

An  array  of  two  gages  was  placed  in  front  of  the  wall  at  the  20  foot  location  in  the  wave 
tank  and  another  pair  of  gages  was  placed  behind  the  wall  at  the  80  foot  section  of  the 
tank.  Both  pairs  were  spaced  two  feet  apart  for  the  0.5  Hertz  wave  tests  and  one  foot 
apart  for  the  0.7  Hertz  wave  tests.  The  pair  of  gages  closest  to  the  wave  maker  was  used 
for  an  array  from  which  both  incident  and  reflected  wave  spectra  could  be  resolved  while 
the  rear  pair  of  gages  permitted  the  resolution  of  transmitted  spectra  as  well  as  the 
spectra  of  waves  reflected  from  the  beach  at  the  end  of  the  wave  tank.  Capacitance 
gages  were  again  used  to  measure  front  and  rear  water  levels  at  the  wall.  The  wall  was 
subjected  to  the  same  set  of  wave  spectra  at  both  wall  penetrations,  requiring  the 


generation  of  only  two  complete  spectra:  one  with  a  modal  frequency  of  0.5  Hertz  and 
the  other  at  a  modal  frequency  of  0.7  Hertz. 
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EXPERIMENTAL  RESULTS 

Preliminary  testing  was  conducted  to  evaluate  the  wall  design  and  functionality  of 
the  support  frame.  Tests  that  were  run  yielded  valuable  ideas  for  design  improvements 
and  indicated  that  quality  data  could  be  obtained  with  some  alterations  to  the  original 
model  design.  Full  time  testing  was  then  initiated  and  completed  and  the  results  are 
described  below. 


WAVE  TRANSMISSION  -  REGULAR  WAVES 

The  experimental  setup  for  this  portion  of  testing  was  the  same  as  previously 
described  and  as  is  shown  in  Figure  16.  The  model  was  subjected  to  waves  of  frequency 
0.5,  0.7,  0.9,  and  1.1  Hertz  and  usually  three  different  steepnesses  of  1/15,  1/20,  and  1/40 
at  each  wave  height.  Waves  with  of  0.5  or  0.7  Hertz  were  not  examined  at  the  1/15 
steepness  for  wall  penetrations  of  2.5  and  3  feet  because  laboratory  facilities  were  not 
capable  of  handling  waves  of  that  magnitude.  Wave  transmission  was  determined  by 
recording  incident  wave  height  with  gage  #1  and  transmitted  wave  height  with  gage  #4 
(Figure  16).  For  both  incident  and  transmitted  wave  records,  five  complete  wave  cycles 
were  identified  and  the  five  wave  crests  and  troughs  were  averaged  to  determine  the 
experimental  wave  height. 

Wave  transmission  was  measured  for  all  testing  sequences  and  then  compared 
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with  the  modified  theory  for  wave  transmission.  The  results  are  shown  below  in  Figure 
21.  This  figure  compares  experimental  values  ofK,  to  theoretical  values  computed  using 
the  modified  theory  for  K, .  In  Figure  21,  data  points  would  fall  along  a  straight  line 
having  a  1:1  slope  if  experimental  results  perfectly  matched  theory.  Figure  21  indicates 
that  agreement  between  experimental  data  and  theoretical  prediction  is  generally  good 
with  some  variations  in  results  for  different  wave  frequencies. 

As  explained  previously,  waves  with  smaller  frequencies  have  longer  wave 
lengths.  As  a  result  these  waves  have  a  deeper  wave  base  and  produce  more 
transmission.  Under  testing  conditions,  a  0.5  Hertz  wave  had  a  wave  length  of  almost  20 
feet  and  thus  a  wave  base  of  10  feet.  This  was  a  large  enough  value  to  be  strongly 
influenced  by  changing  the  penetration  of  the  wall.  Note  in  Figure  21  the  three  distinct 
groups  of  transmission  values  corresponding  to  the  three  tested  wall  penetrations.  The 
smallest  values  of  transmission  for  the  0.5  Hertz  waves  (K,  of  approximately  0.3  to  0.4) 
correspond  to  the  greatest  values  of  wall  penetration  (w  =  3  feet).  At  0.5  Hertz,  theory 
tended  to  underestimate  K, .  The  error  between  data  and  theory,  defined  as 


Kt  exp  -Ktpred 
Ktpred 


(47) 


ranged  from  5  to  21  %  for  0.5  Hertz  waves  using  Wiegel  theory  or  4  to  23  %  using  modified 
theory.  The  0.7  Hertz  waves  had  a  wavelength  of  about  10  feet  and  a  wave  base  of 
approximately  5  feet.  These  waves  produced  some  transmission  with  a  K,  between  0.1  and 


Kt  Experimental 
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Comparison  of  Modified  Theory  and 
Experimental  Data  for  Kt 

Line  of  Perfect 


Figure  21.  Comparison  of  modified  theory  predictions  and  experimental  values  of  the 
coefficient  of  transmission,  K,. 
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0.2.  Theory  tended  to  slightly  overestimate  transmission  for  0.7  Hertz  waves  but  was  more 
accurate  overall.  Modified  theory  had  an  error  of  1  to  22%  while  Wiegel  error  ranged  from 
42  to  58%. 

Little  useful  data  was  obtained  from  the  0.9  and  1.1  Hertz  waves.  The  0.9  Hertz 
waves  were  approximately  6  feet  long  and  thus  were  nearly  completely  reflected  by  the 
wave  barrier.  Scatter  in  the  data  is  greatest  at  this  frequency  because  predicted  and 
actual  values  of  wave  transmission  were  so  small.  Additionally,  random  noise  and 
inherent  error  in  the  system  prevented  accurate  measuring  of  the  extremely  small 
transmitted  wave  heights.  The  cluster  of  1 . 1  Hz  points  at  the  origin  is  the  result  of  these 
waves  being  nearly  perfectly  reflected  by  the  test  model  at  all  tested  penetrations.  These 
data  series  served  mainly  as  a  check  in  wave  interaction  with  an  effectively  “full  depth” 
wall  since  little  or  no  energy  from  these  waves  passed  beneath  the  wave  barrier. 

The  results  of  comparing  experimental  data  from  this  project  with  prediction 
from  Wiegel  (1960)  theory  is  shown  in  Figure  22.  For  most  test  conditions,  Wiegel’s 
theory  predicts  much  larger  coefficients  of  transmission  than  are  actually  observed. 

These  results  agree  with  an  observation  by  Mattson  and  Cederwall  (1976)  which  says, 
“...Wiegel’s  equation  tends  to  give  too  large  values  of  K,  in  the  interval  of  interest  where 
K,  <  0.6.”  When  considered  over  the  tested  range,  the  difference  between  these  two 
methods  is  significant,  ranging  from  15  to  greater  than  200  %.  As  a  result  of  the 
comparisons  between  Figures  21  and  22,  it  may  be  concluded  that  the  modified  theory 
provides  a  much  better  prediction  of  wave  transmission  than  the  Wiegel  (1960)  theory. 
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Comparison  of  Experimental  Data  and 
Wiegel's  Theory  for  Kt 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7 

Kt  from  Wiegel’s  Theory 


Figure  22.  A  comparison  of  Wiegel’s  theoretical  predictions  and  experimental  values  of 
the  coefficient  of  transmission,  IQ. 
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WAVE  .TRANSMISSION  -  IRREGULAR  WAVES 

Experimental  setup  for  random  wave  testing  is  shown  in  Figure  23.  The  pairs  of 
gages  at  20  feet  and  80  feet  were  separated  by  2  feet  for  waves  at  0.5  Hertz  peak 
frequencies  and  were  separated  by  1  foot  for  0.7  Hertz  waves,  as  explained  earlier. 
Testing  was  performed  at  modal  (peak)  frequencies  of  0.5  and  0.7  Hertz  and  at 
penetrations  of  2  and  2.5  feet.  Five  60  second  test  runs  were  made  at  each  combination 
of  penetration  and  peak  frequency  for  a  total  of  20  data  runs  involving  random  waves. 
Data  were  taken  from  six  wave  gages  and  two  force  gages.  The  wave  gages  were 
grouped  in  pairs  with  a  pair  at  20  feet,  a  pair  at  the  wall  (one  inch  in  front  of  and  two 
inches  behind),  and  the  third  pair  at  80  feet.  Both  gages  at  the  wall  and  one  of  each 
extreme  pair  of  gages  were  capacitance  gages  while  the  remaining  two  gages,  one  of  each 
extreme  pair,  were  resistance-type  gages. 

Wave  heights  could  not  be  measured  in  the  same  way  for  irregular  wave  tests  as 
in  regular  wave  tests.  An  irregular  wave  pattern  is  the  result  of  the  random  superposition 
of  sine  curves  at  various  frequencies  and  amplitudes.  Waves  of  different  frequencies 
move  at  different  velocities,  thus  an  irregular  wave  pattern  is  the  result  of  superposition 
of  many  different  waves.  This  superposition  prevents  the  direct  measurement  of  an 
individual  wave’s  height  and  thus  requires  an  alternative  analysis  technique.  In  addition, 
reflected  waves  also  travel  in  a  different  direction  and  further  complicate  wave  height 
measurement.  To  overcome  these  difficulties,  a  method  developed  by  Goda  and  Suzuki 
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(1976)  was  utilized  in  this  project  for  irregular  wave  analysis.  This  method  involves  the 
numerical  analysis  of  signals  obtained  from  a  pair  of  wave  gages  in  close  proximity  to 
each  other.  The  raw  data  is  first  analyzed  using  a  Fourier  Transform.  This  transform  will 
calculate  the  Fourier  coefficients  A  and  B  at  each  gage  for  each  frequency  in  the  wave 
spectrum.  These  coefficients  contain  information  regarding  the  amplitude  and  phase  of 
waves  in  the  spectrum  and,  when  known  at  two  locations,  they  can  be  used  to  separate 
the  incident  and  reflected  waves.  Goda  and  Suzuki  (1976)  developed  a  method  that 
allows  the  prediction  of  incident  and  reflected  wave  amplitudes  at  each  of  the  wave 
frequencies  from  the  Fourier  coefficients  (Figure  24).  After  the  incident  and  reflected 
coefficients  have  been  determined,  further  numerical  analysis  can  then  be  used  to 
calculate  incident  and  reflected  significant  wave  heights  (average  of  1/3  largest  waves  or 
4  times  the  standard  deviation  of  the  water  surface),  which  are  in  turn  used  to  calculate 
the  coefficient  of  reflection,  K, . 

After  the  random  wave  data  was  analyzed  using  the  above  procedure,  average 
values  of  experimental  wave  transmission  (averaged  for  all  five  60  second  test  runs)  were 
compared  to  theoretical  values  predicted  using  the  modal  (peak)  frequency  of  the 
irregular  wave  spectrum.  The  results  of  this  comparison  are  shown  in  Figure  25.  The  0.5 
Hertz  waves  show  fairly  good  correlation  with  theory.  They  have  the  same  slope  as  the 
line  of  theoretical  prediction  but  fall  below  this  line.  This  is  because  these  results  are  not 
for  a  pure  wave  of  a  single  frequency  but  rather  a  wave  spectrum  that  consists  of 
frequencies  ranging  from  approximately  0.4  to  1.0  Hertz.  The  0.5  Hertz  waves  in  the 
spectrum  may  be  reasonably  expected  to  follow  theoretical  predictions,  but  the 
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Further,  we  suppose  that  the  surface  elevations  are  recorded  at  two 
adjacent  stations  of  xj  and  x2  «  xj  +  4£.  The  observed  profiles  of  composite 
waves  will  be 


ni 

n2 

where, 


(n, 


(nn 


+  V 

+  nR) 


X«=Xj 

xcx2 


Aj  cos  at  +  B i sin o  t , 
A2cos  c  t  +  B2sin  o  t , 
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Equation  3  can  be  solved  to  yield  the  estimate  of 
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In  the  calculation,  the  dispersion  relation  of  the  following  is  presumed 
to  hold: 


o2  «  gk  tanh  kh. 


(6) 


Figure  24.  Goda  and  Suzuki  (1976)  method  of  irregular  wave  analysis. 


Comparison  of  Modified  Theory  and 
Experimental  Data  for  Random  Waves 


Kt  from  Modified  Theory 


a  0.5  Hz  modal ,  w  =  2.5'  •  0.7  Hz  modal,  w  =  2.5' 
■  0.5  Hz  modal,  w  =  2' 


Figure  25.  Comparison  of  K,  Experimental  with  K,  Theoretical  for  irregular  wave  tests 
on  wave  barrier. 
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components  of  the  spectrum  in  the  higher  frequencies  will  have  much  smaller  amounts  of 
transmission.  Thus,  the  overall  coefficient  of  transmission  for  the  wave  spectrum  will  be 
smaller  than  modified  theory  suggests  for  regular  waves  for  a  0.5  Hertz  wave.  While  this 
difference  is  pronounced,  it  is  not  large  due  to  the  fact  that  only  a  small  amount  of  the 
spectral  energy  will  be  distributed  to  those  waves  of  frequency  farthest  from  the  peak. 
This  reasoning  explains  the  0.5  Hertz  test  data,  but  at  this  time  the  0.7  Hertz  data  trends 
cannot  be  accounted  for. 


WAVE  FORCES  -  REGULAR  WAVES 

Wave  forces  due  to  regular  waves  were  recorded  at  the  same  time  as  data  for 
transmission  of  regular  waves  and  thus  the  setup  and  testing  procedures  follow  those  as 
outlined  previously  for  the  wave  transmission.  Forces  were  recorded  individually  near 
the  top  and  bottom  of  the  wave  barrier  using  two  force  gages  (Figure  18)  and  these 
measurements  were  then  summed  to  determine  a  single,  total  force  on  the  center  test 
section.  These  total  forces  were  then  compared  to  predictions  from  the  modified  theory 
for  forces  (Equation  41)  and  the  results  are  shown  in  Figure  26.  Forces  in  the  downtank 
direction  (resulting  from  a  wave  crest  at  the  wall)  are  positive  and  are  denoted  by  a 
triangle.  Forces  obtained  when  a  trough  is  at  the  wall  are  negative  and  denoted  by  an 
oval.  The  absolute  values  of  the  experimental  forces  were  used  in  this  graph.  From 
Figure  26,  it  is  apparent  that  the  modified  theory  for  computing  forces  is  generally 
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Evaluation  of  Modified  Theory  with 
Phase  Shift  and  Second  Order  Force 


0  10  20  30  40  50  60  70  80  90 

Theoretical  Force  (lbs) 


Figure  26.  Comparison  of  experimental  data  to  modified  force  theory  predictions, 
waves. 
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accurate.  Overall,  the  new  theory  provides  relatively  small  errors  in  force  prediction. 

This  is  especially  true  for  the  largest  forces  measured  (from  30  to  80  lb  due  to  0.7  and  0.5 
Hertz  waves)  where  errors  were  less  than  20%. 

On  average,  the  theory  is  slightly  conservative.  By  applying  linear  regression 
theory  to  the  experimental  data,  it  was  determined  that  experimental  force  was  98.7%  of 
the  experimental  force.  Most  of  the  largest  error  occurred  for  0.9  and  1 . 1  Hertz  waves 
where  measured  forces  were  fairly  small. 

The  testing  matrix  for  regular  wave  tests  developed  two  distinct  situations  when 
predicting  forces.  The  first  situation  applies  to  0.5  and  0.7  Hertz  waves.  These  waves 
are  long  enough  to  transmit  some  of  their  energy  under  the  barrier  and  the  coefficient  of 
transmission  for  these  waves  has  an  influence  on  how  much  force  they  exert  on  the 
barrier.  These  waves  generally  produced  forces  in  excess  of  25  pounds  on  the  barrier 
and  the  modified  theory  had  to  account  for  the  effects  of  the  transmission  and  reflection 
of  these  waves.  On  the  other  hand,  the  transmission  of  the  0.9  and  1 . 1  Hertz  waves  was 
almost  always  equal  to  zero.  The  predicted  coefficients  of  transmission  and  reflection 
therefore  played  no  role  in  predicting  the  forces  on  the  barrier  at  these  frequencies.  The 
predictions  of  the  modified  equation  were  still  accurate,  however.  This  shows  the  wide 
range  of  applicability  of  the  modified  theory  because  its  accuracy  does  not  vary  with 
different  ratios  of  wave  barrier  penetration  relative  to  incident  wave  length  or  any  other 
factors. 

Figure  27  shows  the  results  of  using  the  NAVFAC  method  (NAVFAC,  1982)  to 
compute  the  force  on  the  model  under  test  conditions.  The  NAVFAC  force 
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Comparison  of  Experimental  Forces 
and  NAVFAC  Theory  for  Forces 


NAVFAC  Prediction  (lbs) 


Figure  26.  Comparison  of  NAVFAC  predicted  forces  to  experimental  data  for  regular 
waves. 
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predictions  differ  considerably  from  those  of  the  modified  theory.  In  this  case,  forces  are 
overstated  by  a  factor  of  two  to  six  (errors  of 200%  to  600%).  This  large  error  is  a 
result  of  incorrect  assumptions  used  in  force  computation  by  the  NAVFAC  manual.  The 
largest  source  of  error  is  due  to  the  assumption  that  when  a  crest  is  acting  on  one  side  of 
a  wall,  a  trough  will  be  acting  on  the  other  side.  For  waves  with  significant  amounts  of 
transmission  such  as  the  0.5  and  0.7  Hertz  waves  used  in  this  project,  this  assumption  will 
result  in  the  overestimation  of  force  by  a  factor  of  two  to  three.  This  error  is  due  to  the 
large  pressure  difference  that  results  with  an  elevated  water  surface  at  the  front  of  the 
wall  and  a  depressed  water  surface  behind  the  wall  (assuming  a  crest  at  the  front).  As 
noted,  observations  of  water  levels  recorded  in  this  study  show  that  there  is  always  a  90 
degree  phase  shift  in  water  levels  from  the  front  to  the  back  of  the  wall.  When  a  crest  is 
at  the  front,  the  zero  crossing  of  the  wave  is  at  the  rear,  i.e.  the  water  surface  at  the  rear 
is  at  the  still  water  level. 

While  this  incorrect  assumption  can  be  the  source  of  a  large  amount  of  error, 
another  incorrect  assumption  can  generate  even  larger  factors  of  error.  The  NAVFAC 
method  approximates  dynamic  pressure  due  to  incident  waves  as  a  straight  line  when  in 
actuality  it  is  a  hyperbolic  relationship  according  to  the  linear  wave  theory  for  dynamic 
pressure  (Equation  15).  The  NAVFAC  method  then  further  compounds  this  error  for 
deep  water  waves  (such  as  the  0.9  and  1.1  Hertz  waves  used  in  this  experiment)  by 
assuming  that  the  dynamic  pressure  is  equal  to  zero  at  the  sea  floor.  In  reality,  dynamic 
pressure  is  equal  to  zero  at  the  wave  base.  For  waves  with  the  wave  base  at  a  significant 
height  above  the  sea  floor,  this  assumption  can  result  in  large  errors  when  the  pressure  is 
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integrated  over  the  entire  depth. 


UNDERWALL  FLOW  TESTING 

Fluid  velocities  in  front  of  and  under  the  wall  were  measured  and  recorded  using 
an  Acoustic  Doppler  Velocimeter  at  the  positions  shown  in  Figure  20.  This  series  of 
testing  utilized  regular  waves  having  a  frequency  of  0.5  Hertz,  a  wave  steepness  of  1/40, 
and  a  wall  penetration  of  2.5  feet.  Particle  velocities  under  the  wall  were  examined  to 
gain  a  better  idea  of  the  velocity  and  direction  of  flow  underneath  the  wall. 

Over  a  series  of  five  waves,  horizontal  and  vertical  particle  velocities  were 
recorded.  The  magnitude  of  the  velocity  was  then  calculated  for  each  pair  of  velocities 
and  over  a  series  of  five  waves  the  recorded  velocities  and  magnitudes  were  averaged. 
Additionally,  a  least-squares-fit  analysis  was  performed  to  determine  the  average  error  of 
each  data  run.  A  least-squares-fit  error  analysis  examines  an  experimental  data  series 
such  as  that  of  particle  velocity  versus  time.  Through  numerical  analysis,  amplitude 
coefficients  are  then  calculated  which,  when  multiplied  by  sine  and  cosine  components, 
will  approximate  an  “ideal”  curve  through  a  data  series.  The  error  at  any  point  can  be 
calculated  by  comparing  the  value  of  this  “ideal”  curve  at  a  point  with  the  value  of  the 
experimental  data  at  the  same  point.  This  process  was  repeated  at  each  grid  point  by 
generating  another  set  of  waves  after  moving  the  ADV  probe  to  a  location  shown  on  the 
grid. 
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Particle  velocities  were  plotted  as  a  function  of  time  and  as  a  function  of  depth  in 
order  to  compare  flow  at  different  points.  As  shown  in  Figure  28,  water  particle  velocity 
under  the  wall  did  depend  strongly  on  depth.  There  is  a  significant  departure  between 
particle  velocity  as  predicted  by  theory  and  actual  particle  velocity  in  the  area 
immediately  below  the  wall.  This  departure  exists  because  particles  at  this  area  cease  to 
travel  in  their  usual  elliptical  orbits.  There  is  a  net  surge  of  water  back  and  forth  under 
the  wave  barrier  due  to  the  large  pressure  difference  between  the  front  and  back  of  the 
wall.  Water  will  flow  from  an  area  of  higher  pressure  to  lower  pressure.  Thus,  when  a 
wave  crest  is  at  the  front  of  the  wall  dynamic  pressures  are  at  their  maximum,  water 
begins  to  flow  from  the  front  of  the  wall  to  the  back  of  the  wall.  When  there  is  a  trough 
at  the  front  of  the  wall  the  lower-pressure  region  is  now  in  front  and  water  will  flow  back 
to  the  front  side  of  the  wall.  Figure  29  shows  the  turbulence  in  this  back-and-forth  flow 
beneath  the  bottom  of  the  wall.  Figure  30  is  a  graph  of  particle  velocities  at  a  greater 
depth  than  Figure  29  and  shows  the  general  smoothness  of  flow  at  greater  depths.  Flow 
at  shallower  penetrations  is  much  more  turbulent.  This  can  be  shown  in  two  ways.  First, 
the  regular  horizontal  and  vertical  particle  flow  is  more  irregular  at  shallower  depths. 
Second,  there  is  a  significant  amount  of  motion  in  the  Z- ,  or  cross-tank,  direction  at 
shallower  depths.  Particle  motions  should  be  in  only  two-dimensions  according  to  linear 
wave  theory.  The  presence  of  a  vortex  as  a  result  of  turbulent  flow  could  induce  particle 
velocities  in  the  Z-direction,  however.  It  is  this  large  amount  of  cross-tank  flow  that  also 
helped  to  determine  the  size  of  the  induced  vortex  in  the  next  part  of  the  project. 


Horizontal  Particle  Velocity  (ft/s) 
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Decay  of  Experimental  and  Theoretical 
Particle  Velocity  (U)  with  Depth 


—  Theoretical  U  -®-  Experimental  U 


Figure  28.  Graph  of  experimental  and  theoretical  horizontal  velocity  (with  wave  barrier 
annotated)  versus  depth. 
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Underwall  Velocity  Profile 
3"  Below  Bottom  of  Wall 


Time  (s) 

—  Ux  — —  Vy  —  Z  - H 


Figure  29.  Graph  of  underwall  particle  velocity  versus  time  at  a  penetration  of  33  inches. 


Underwall  Velocity  Profile 
19"  Below  Bottom  of  Wall 


Figure  30.  Graph  of  underwall  particle  velocity  versus  time  at  a  penetration  of  49  inches. 


Water  Level  (in) 
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VORTEX  FLOW  TESTING 

Vortex  flow  tests  were  performed  under  the  same  conditions  as  the  underwall 
flow  tests.  The  goal  of  this  part  of  testing  was  to  map  the  extent  and  magnitude  of  the 
vortex  that  forms  from  the  irregular  flow  under  the  wall.  Data  from  this  phase  were 
analyzed  with  the  same  methods  as  used  in  underwall  flow  analysis.  A  graph  of  least- 
squares-fit  error  at  each  testing  point  is  shown  in  Figure  3 1 . 

This  error  was  calculated  by  first  computing  the  least-square  fit  of  a  0.5  Hertz 
sine/cosine  waveform  to  the  velocity  signals  recorded  during  testing.  The  error  between 
each  measured  data  point  and  the  best-fit  sine/cosine  curve  was  then  calculated  for  all 
points  and  averaged.  This  mean  error  is  the  number  next  to  each  measurement  point  in 
Figure  31. 

The  calculation  of  the  mean  error  can  be  used  to  map  the  vortex  if  it  assumed  that 
the  deviation  from  the  best-fit  sine/cosine  function  is  the  result  of  turbulent  flow.  The 
presence  of  a  baseline  error  seems  to  validate  this  assumption.  If  the  points  45"  and 
greater  below  the  wall  are  examined,  there  appears  to  be  an  average  error  of 
approximately  0.01.  By  assuming  this  error  to  be  the  “background”  or  “baseline”  error,  a 
reference  is  gained  for  gauging  the  turbulence  of  the  vortex  in  different  areas. 

The  error  in  Figure  31  is  largest  in  the  area  surrounding  the  bottom  of  the  wall. 
The  average  error  then  decreases  as  points  farther  from  this  “pocket”  are  considered. 

This  corresponds  with  fluid  dynamics  theory  because  vortices  are  supposed  to  lose 
energy  as  they  move  away  from  their  point  of  origin.  Thus  the  area  where  the  vortices 
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are  the  strongest  and  would  therefore  generate  the  most  error  is  immediately  surrounding 
their  point  of  origin,  which  in  this  case  is  the  bottom  of  the  wall.  As  points  farther  from 
this  origin  are  examined,  the  error  decreases,  as  would  the  energy  of  any  turbulent  flow 
as  it  approaches  those  points. 


Depth  (in) 
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CONCLUSIONS 

From  the  results  of  this  study  on  vertical  wave  barriers,  the  following  conclusions 
may  be  drawn: 

-  It  is  necessary  to  account  for  wave  reflection  when  attempting  to  calculate  wave 
transmission  and  forces.  Wiegel  theory  does  not  account  for  reflection  and  its  effects  on 
particle  velocity  and  wave  pressures.  This  causes  Wiegel  theory  to  overestimate  wave 
transmission. 

-  The  modified  theory  for  transmission  is  a  valid  improvement  on  existing  methods  of 
calculating  wave  transmission.  The  experimental  data  corresponded  well  with  predicted 
values  of  transmission. 

-  There  is  a  90-degree  phase  shift  in  the  water  level  between  the  front  and  back  of  a 
vertical  wave  barrier.  This  phase  shift  ensures  that  when  there  is  either  a  crest  or  a 
trough  at  the  front  of  a  wave  barrier,  the  water  behind  the  wall  will  be  at  the  still  water 
level. 

-  First  and  second  order  dynamic  forces  and  the  90-degree  phase  shift  must  be  accounted 
for  in  the  devlopment  of  a  theory  to  predict  forces  on  a  wave  barrier.  The  modified 
theory  for  wave  forces  includes  these  considerations  and  is  a  valid  theory  that  can  be 
expected  to  predict  wave  forces  on  a  vertical  wave  barrier  accurately. 

-  The  existing  NAVFAC  method  for  computing  forces  on  wave  barriers  is  extremely 
conservative.  The  assumption  of  an  1 80-degree  phase  shift  is  incorrect,  as  is  the 
assumption  that  dynamic  pressures  will  extend  to  the  seafloor  instead  of  the  wave  base. 
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-  The  back-and-forth  flow  of  water  beneath  a  wave  barrier  is  turbulent  and  induces  the 
formation  of  vortices.  These  vortices  act  primarily  in  an  area  slightly  above  the  bottom  of 
a  wave  barrier. 

-  When  analyzing  vertical  wave  barriers,  an  assumption  of  conservation  of  mass  should  be 
used  instead  of  conservation  of  energy.  The  formation  of  vortices  from  turbulent  flow 
confirms  that  energy  is  not  conserved  during  water  flow  beneath  a  wave  barrer. 

-The  modified  theories  for  transmission  and  wave  forces  significantly  improve  the 
prediction  of  wave  transmission  and  the  forces  on  a  wave  barrier.  In  addition,  the 
modified  theory  for  wave  forces  is  much  simpler  and  easier  to  use  than  the  NAVFAC 
method  for  predicting  force. 
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